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• The atmospheric neutrino-fiux is calculated over the wide energy range from 1 GeV 

' to 3,000 GeV for the study of neutrino-physics using the data from underground neutrino- 

\^ , detectors. The uncertainty of atmospheric neutrino-fluxes is also discussed. A brief comment 

■ is made to interpret the anomaly in terms of neutrino oscillations. 

m ■ 

> : 

ff^ , §1. Introduction 

(N ■ 

. Since the observation of atmospheric neutrinos by the Kamiokande group B, 

T-H ', some of the underground detectors a W have confirmed that the flux ratio of neutrino 

species (f^ + f^) / (fg + fe) is significantly different from the expected value, although 
. the situation is still controversial a>a>. The atmospheric neutrino anomaly may 

have a crucial importance in particle physics, sine it can be interpreted in terms of 
neutrino oscillaticins with a large mixing angle and a typical mass squared difference 
of O (10^^ eV^) ePQ) eI^b). The observation of the zenith angle variaton of the double 
(— I I ratio (/i/e)rfato/(/i/e)MC at multi-GeV energies is also suggestive li^P. It is therefore 

' important to calculate atmospheric neutrino fluxes precisely. 

•'— j . Atmospheric neutrino fluxes have been calculated from the incident beam rvf 

r> I primary cosmic rays by Volkovatil', Mitsui et alllH^, Butkevich et alE^), and LiparilliP 

' mainly for high energies (from around 1 GeV to above 1 00 00 GeV). Gaisser et 

alEJ, Barr et alt§, Bu gaev and NaumovlllP, Lee ans KohGJ*, and Honda et alllHP 
made a detailed calculation of the atmospheric neutrinixfluxes for low energies from 
the primary cosmic rays. On the other hand, Perkins c^) calculated the low energy 
atmospheric neutrinos using ^-flux observed at high altitude. 

In this paper, we report the detailed calculation of the atmospheric neutrino 
fluxes in the energy range from 30 MeV up to 3,000 GeV, corresponding to the ob- 
servation range of underground neutrino detectors til-*. We also discuss the possiblity 
to interpret the anomaly in terms of neutrino oscillations. 



*' Talk presented by S. Midorikawa, to appaer in the Proceedings of Yukawa International 
Seminar '95: From the Standard Model to Grand Unified Theories. 
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§2. Primary Cosmic Ray Fluxes 
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Primary cosmic ray fluxes are relatively well known in the low energy region(< 
lOOGeV), by which the low energy atmospheric ly-Quxes (< 3 GeV)are mainly pro- 
duced. Webber and LezniakBHP have compiled the energy spectrum of the cosmic 
rays for the hydrogen, helium, and CNO nuclei in the range 10 MeV ~ 1,000 GeV 
for three levels of solar activity A similar compilation has been made by others for 
hydrogen and helium nuclei, which agrees well with that of Webber and Lezniak. 

The geomagnetic field 
determines the minimum en- 
egy with which a cosmic ray 
can arrive at the earth. For 
the cosmic ray nucleus, the 
minimum energy of cosmic 
rays arriving at the earth 
is determined by the mini- 
mum rigidity (rigidity cutoff) 
rather than the minimum 
momentum. The value of the 
rigidity cutoff for the actual 
geomagnetic field can be ob- 
tained from a computer sim- 
ulation of cosmic ray tra- 
jectories. If a cosmic ray 
particle can reach the earth, 
the antiparticle with an op- 
posite momentum can escape 
from the earth. We launch 
antiprotons from the earth, 
varying the position and di- 
rection. When a test parti- 
cle with a given momentum 
reaches a distance of 10 times 
of the earth's radius, it is 
assumed that the test parti- 
cle has escaped from the geo- 
magnetic field. The rigidity cutoff at Kamioka site is shown as the contour map in 
Fig. 2. 

Cosmic rays with energy greater than 100 GeV, which are responsible for ^10 GeV 
atmospheric neutrino fluxes, are not affected by solar activity and by geomag- 
netic field. There are few measurements of the cosmic ray chemical composition 
at these energies, especially above 1 TeV. We compiled the available data and and 
parametrized the observed fluxes for > 100 GeV with a single power function, and 
show the result in table |. We treated bound nucleons at these energies as indepen- 
dent particles, and estimated the primary nucleon spectrum. 
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Fig. 1. Observed fluxes of cosmic ray protons, helium nu- 
clei, and CNOs from the compilation of Webber and 
Leziak^-*. Solid lines are our parametrization for solar 
mid, dashed lines for solar min., and dotted lines for 
solar max. 
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Fig. 2. The contour map of cutoff-rigidity for the u arrival directions at Kamioka. Azimuthal angles 
of 0°, 90°, 180°, and 270° show directions of south, easth, north, and west respectively. 

Table I. Compiled cosmic ray spectrum in the form: A{E/100 GeV)''. 



Nucleus 


A 


7 


H 


(6.65 ±0.13) X 10"^ 


-2.75 ± 0.020 


He 


(3.28 ±0.05) X 10"^ 


-2.64 ±0.014 


CNO 


(1.40 ±0.07) X 10"" 


-2.50 ±0.06 


Ne-S 


(3.91 ± 0.03) X lO"'"' 


-2.49 ±0.04 


Fe 


(1.27 ±0.11) X 10"^ 


-2.56 ±0.04 



§3. Production and Decay of hadrons 

As cosmic rays propagate in the atmosphere, they produce vr's and K's in in- 
teractions with air nuclei. These mesons create atmospheric i^'s when they decay as 
follows: 



A-cr + ^air * i K , K , ■ ■ ■ 

7r+ ^ + 

/i+ ^ e+ + l^e + t'/i 
vr" fi~ + (3-1) 

11' ^ e' +Ve + I^fj. 



The calculations of the cosmic ray protons with air nuclei consists of a number 
of Monte Ilarlo codes corresponding to different primary energies. We employed the 
NUCRINE3) Monte Carlo code for the hadronic interaction of cosmic rays for Ei^ < 
5 GeV, and LUND code - FRITIOF version 1.60 and JETSET version 6.3© - 
for 5 GeV < Eiab < 500 GeV. Above 500 GeV, the original code developed by 
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Kasahara et al (COSMOS)© was used. The K/tt ratio is taken 7 % at 10 GeV, 
11 % at 100 GeV, and 14 % at 1,000 GeV in laboratory energy. 

We consider all the decay modes of vr and K mesons but for rare ones. We have 
ignored charmed meson production, since the contribution of charmed particle to 
atmospheric neutrinos becomes sizable only for Ei^>WO TeV. 

In the two body decay of charged tt's and iC's, the resulting fi^ is fully polarized 
against (toward) the direction of fi motion in the charged vr or rest frame. We 
took into account the muon polarization effect in the subsequent decay following 
HayakawacIP. We applied the discussion in Ref. 28 for the polarization of //'s from 
the K^ii, decay. The small angle scattering of /x's in the atmosphere reduces the /x 
polarization. This depolarization effect was also evaluated in Ref. 27 as of the order 
of 21 MeY/vp, where v and p are velocity and momentum of ^u's respectively. 

§4. Atmospheric Neutrinos 

At low energies, the rigidity cutoff has a significant directional variation. In 
the one-dimensinoal approximation which we adopted, we expect larger z^-fluxes 
from the low rigidity cutoff directions and a smaller z^-fluxes from the high rigidity 
cutoff direction. In the actual case, however, it may be difficult to observe these 
variations. There is a smearing effect of direction in the i/-detector. When a low 
energy u ( < 3 GeV) creates a charged lepton by a quasi-elastic process, the lepton has 
a typical angle of 50 - 60° from the u direction. Thus the directional dependence of 
atmospheric neutrino flux is small for lower enregy neutrinos, especially when they 
are observed in the detector. We present in Fig. 3 the atmospheric neutrino fluxes 
averaging over all directions together with other calculations. 

In Fig. 4, we show the flux ratio by i/-species along with those of other authors. 
Although the calculation method and some of physical assumptions are different 
among these authors, the ratio {ue + ^'e)/{^fi + i^fi) is very similar each other. The 
relatively large difference in v'e/'^e among them may reflect the difference of calcula- 
tion scheme and/or the physical assumptions. 

§5. Uncertainties 

The systematic error in the calculation of atmospheric j^~fluxes comes mainly 
from the incompleteness of the knowledge of the primary cosmic ray fluxes. Even 
at low energies, where the primary cosmic ray fluxes are rather well studied, it is 
difficult to determine the absolute value due to the uncertainties in the instrumental 
efficiency (~ 12%) and exposure factor (2 - 3%). In our compilation, the error in the 
fit is ~ 10 % for the nucleon flux at 100 GeV and ~ 20 % at 100 TeV. Assuming 
~ 10% uncertainty below 100 GeV, the systematic error in the atmospheric i^-fluxes 
is estimated to be ~ 10% at < 3 GeV, increasing to ~ 20 % at 100 GeV, and remains 
almost constant up to 1,000 GeV. 

The interaction model is another source of systematic errors. In our comparison, 
the agreement of the LUND model and the COSMOS code with the experimental 
data is <, 10 %. The authors of the NUCRIN code claim that the agreement is within 
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Fig. 3. The atmospheric i^-fluxes multiplied by for the Kamioka site at solar mid. (solid line). 
BGS are from Ref. [16], BN from Ref. [17], and LK from Ref. [18]. The dotted line is the resuh 
from the calculation for high energy without the rigidity cutoff, and averaged over all directions. 
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Fig. 4. The flux ratio of i/-species calculated for Kamioka. BGS are from Ref. [16], BN from 
Ref. [17], and LK from Ref. [18] as before. 



10 - 20 The hadronic interaction below 5 GeV contributes at most 5 % to the 

production of atmospheric z^-fluxes at 1 GeV. The systematic error caused by the 
hadronic interaction model is estimated to be ~ 10 % above 1 GeV. 

One-dimensional approximation which we have adopted is justified only at high 
energies. It is expected to be accurate above 3 GeV. Since the calculation of rigidity 
cutoff is very simplified in this scheme, this may result in a systematic error in the 
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absolute value of the atmospheric neutrino fluxes of 10 - 20 % at 100 MeV and 5 % 
at 1 GeV. 

The statistics of the Monte Carlo calculation also causes an error in the atmo- 
spheric neutrino fluxes. The uncertainty due to the statistics is estimated to be ^5 % 
up to 100 - 300 GeV for i/^ and z/^, and up to 30 - 100 GeV for and v'e, depending 
on the zenith angle. The errors increase to ~ 10 % at the highest energy for each 
kind of i^'s. 

Combining all the systematic and non-systematic errors, the total error is es- 
timated as 15 % from 1 GeV to 100 GeV, and 20 - 25 % at the highest energy 
in our calculation. However, the error of the species ratio is smaller than that of 
the absolute value, since the z^~species ratio is not affected much by the uncertainty 
of primary fluxes and the calculation scheme. It is estimated to be <, 10 % below 
100 GeV for h'/i) and ^ 5 % below 30 GeV for {u^ + i^^)/(fe -|- I'e)- These errors also 
increase to 10 - 15 % at the highest energies in our calculation. 

§6. Neutrino Oscillations 

We compare the Kamiokande data0^ with the theotetical calculations. As in 
our previous paper cP, we define. 



< EaaaFfs >=J2 £a{Ea)(ya{Eu,Ea)Fp{E^,, Q^) p{h)dEadEydQydh, (6-1) 

where ea{Ea) is the detection efficiency for an a-type charged lepton with energy Ea, 
Get the differencial cross section of i^q,, Fp{E^, Qy) the incident flux with energy 
El, and zenith angle Qy . 

Insted of the number A'^^^'^ and A^^* of the observed electron and muon events, 
we use the ratios defined as follows: 



< e^a^F^ > ' < EeCTeF^ > 



h = — ^^^rf-, (6-2) 



and 



= K < eJa^F, > = f'Np^ ^'''^ 

N N"^^ 
K<e,a,F,>- NMC^ 4j 

where {k = [number of nucleons] x [time]), N^^^ and Njf^ are the expected 
numbers of electron and muon events from the Monte Carlo calculations respectively. 

If there is no 'atmospheric neutrino anomaly', the data would point to (/2, 1) 
in the (Ue,U^) plain. Note that the effects of flux models, geomagnetic cutoffs, 
and the detection efficiencies are all inluded in the values of /i and /2. We find 
/i ~ /2 ~ 0.473 = / from our neutrino fluxes and the detection efficiency of the 
sub-GeV data at Kamiokande. 
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If there are neutrino oscillations, the ratios Ue and become 



Ue = f2< P{Ve ^Ve)> + < P{v^ ^ V,) > (6-5) 
U^ = < P{U^ ^ l/^) > +/l < Pil^e ^l^^,)>, (6-6) 

where the brackest means the average over the distances and energies of neutrinos: 

1 

< P(l//3 ^ l^a) >= — ^T— X 

f ea{Ea)ac,{E^,Ec,)Fp{E^,Q^)P{up ^ u^)p{h)dEadE^dn^dh. (6-7) 



Table II. the values of Ue, Uf_i, and Ue/U^i calculated from both the sub- and multi-GeV data for 
the Kamioka site with an exposure of 7.7 ktn ■ yr. The value / is taken to be / = 0.473. 



Energy 






Ue/U^ 


Sub-GeV 


0.518 ±0.084 


0.656 ±0.103 


0.790 ± 0.072 


Muhi-GeV 


0.697 ±0.126 


0.832 ±0.144 


0.838 ±0.111 



We summarize in Table II the values of Ue, U^, and Ue/U^ which are obtained 
from both the sub- and multi-GeV data with 7.7 ktn ■ yr. We used only single ring 
events for the analysis of sub-GeV data. For the multi-GeV data, we used both 
single and multi ring events, and evaluated Ue using Eq. 6-4 with /2 = / obtained 
from the sub-GeV data so that we can compare both data directly. We combine the 
statistical and systematic errors of ~ 15% including the uncertaity in the neutrino 
cross section. 

We show in Fig. 5 the regions allowed by the Table II together with the region 
{Ue, Ufj) allowed by three neutrino oscillations assuming /i = /2 = /. 

From this figure, we find that the sub-GeV and muti-GeV data are marginally 
consistent each other, and that it is possible to explain the anomaly in terms of three 
neutrino oscillations. Although it may be premature to deduce what mode dominates 
oscillatios, the multi-GeV data seem to sugest Ve < — > oscillations, while the sub- 
GeV data prefer < — > oscillations. It seems that Ue and U^ grow with energy, 
keeping the C/^/[/e_ratio almost constant. The same conclusion has also been drawn 
by Fogli and Lisic3). This effect may be explained by the improvement of of the 
data and/or calculations. However, if we take both the theory and the experiment 
seriouly, we are lead to a more fascinating conjecture. We can explain both data if 
matter enhanced oscillations would occur between sub- and multi-GeV regions. 

In conclusion, we have shown a detailed calculation of atmospheric neutrinos and 
a new analysis of the anomaly using the recent Kamiokande data. The atmospheric 
neutrino anomaly is eager for further confirmation by another types of neutrino 
expriments such long baseline experiment. 
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Fig. 5. Analysis of the Kamiokandc data. Bold line : Allowed region by three neutrino oscllations. 
Solid line: Kamikande Sub-GeV data. Dashed line: Kamiokande Multi-GeV data. 
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